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Abstract 



Motivated by a recent rather surprising conclusion based on the 1992 PDG 
data on the pion, kaon and lepton decays that if three generations of neutrinos 
are assumed to be massive and mixed, the heaviest neutrino, 1/3, could have 
a mass in the range, 155 MeV < 7713 < 225 MeV, we have analyzed the 
latest 1995 data on the leptonic decays of pion, \x and r with the assumption 
that three generations of neutrinos are massive and mixed. It is shown that 
when the radiative corrections are included and the constraint from partial 
decay widths is imposed, the 1995 data are consistent with three massless 
neutrinos with no mixing. Various limits on the neutrino mass and mixing 
angle implied by the 1995 data are presented together with a critique of the 
previous analysis. 
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I. INTRODUCTION 



In a series of seminal papers M in the early 1980's, Shrock proposed a wide range of 
experimental methods to obtain possible limits on the neutrino mass and associated mixing, 
all based on the precision analysis of the weak interaction data on the decays of the pion, kaon 
and charged leptons, /i and r. In these papers, the comprehensive analysis of decay rates 
and branching ratios of the lepton and meson decays was carried out using the theoretical 
framework with three massive neutrinos and associated mixing. (It is interesting to note 
that as early as in 1961 Bahcall and Curtis @ proposed a similar method based on pion and 
muon decays, even before the discovery of v^Y) At that time, however, available data were 
not accurate enough to provide any significant results on the limits on the neutrino mass 
and mixing angle, in the sense that the limit on the mixing angles were restricted mostly 
for large values of the neutrino masses. In the latter works |||| the limits were further 
improved. 

The best known and often quoted limits on the neutrino mass still come from the analysis 
of spectral shapes in the Kurie plots or other decay kinematics [Q: 

m\ = m(vi) < 5 eV [4] 

m 2 = m(u 2 ) < 270 KeV Q (1) 
m 3 = m{u 3 ) < 24 MeV §. 

The results from these analyses are always presented with the assumption of neutrinos 
with no mixing. A full analysis of the spectral shapes with three massive neutrinos with 
mixing is very much involved and so far no such analysis with satisfactory accuracy has been 
carried out. 

In the meantime, much attention has been focused on entirely different approaches in 
which neutrino mass and associated mixing can be probed indirectly by searching for neu- 
trino oscillation phenomena. The recent activities in this approach include the experimental 
search for reactor and accelerator neutrino oscillations and the study of the solar and at- 
mospheric neutrinos. Although very intriguing indications of massive neutrinos with mixing 
have recently been hinted in the solar and atmospheric neutrino experiments and in the 
LSND experiment, a definitive answer from these experiments is yet to come. 

Recently, Peres, Pleitez and Zukanovich Funchal (PPZ) |/J carried out a comprehensive 
analysis of the existing data on the meson and lepton weak decays assuming that three 
generations of neutrinos are massive and mixed. Their analysis was based on the 1992 
Particle Data Group (PDG) data combined with the latest (in 1993) data on r decays 
||. The results are quite surprising in that the 1992 data on the decay branching ratios would 
be consistent with a finite mass for z/3, i.e. 155 MeV < m.3 < 225 MeV. This mass range is 
significantly larger than the most recent upper limit, m 3 < 24 MeV that was obtained from 
a kinematical analysis of the r decay into five or six pions and v T with no mixing. Moreover, 
PPZ found that the mixing angle (3 which represents mixing between v\ and z/ 3 is also finite 
(11~12 degrees) whereas the mixing angle 7 between v 2 and z/3 was bounded from above, 
thus allowing zero mixing angle. 

Motivated by these rather surprising results, we have carried out a similar analysis of 
the decay rates and branching ratios of the leptonic decays of the pion, /1 and r with the 
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assumption of massive neutrinos with mixing. First we have repeated the PPZ analysis with 
the same set of the data (1992 PDG data) and with the assumptions used by PPZ, confirming 
their results. However, we have also found that their results are significantly modified when 
the constraint coming from the partial decay rates, which PPZ did not use, is imposed. 
The constraint imposed by the decay rate is not an independent one. Instead, it ensures 
that possible fortuitous cancellations in the ratios would not lead to erroneous conclusions. 
Furthermore, we have found that the radiative corrections, which PPZ also ignored, are 
quite important because the accuracy of the data more than warrants the inclusion of the 
radiative corrections in a precision analysis such as this. A similar analysis using the latest 
1995 data shows that the 1995 data is consistent with the picture of three massless neutrinos 
with no mixing. We have been able to set various limits on the neutrino masses and mixing 
angles. 

The plan of our paper is as follows. In Section II, we list the formulas relevant to our 
analysis. (Some details are put in Appendix.) All the data used in our analysis are collected 
in Section III, including the data of 1992 for comparison. In Section IV, we present details 
on our reexamination and critique of the PPZ analysis. In particular, we discuss here what 
would happen to the PPZ conclusions based on the 1992 PDG data, when the radiative 
corrections are included and the constraint from the partial decay rates is imposed. New 
results based on the 1995 data are presented in Section V, and a summary and conclusions 
are given in Section VI. 



II. FORMULAS WITH THREE GENERATION MIXING 



The mixing matrix V in the lepton sector which relates the (weak) interaction eigenfields 
u a (a = e, /i, r) to the mass eigenfields Vi [i = 1, 2, 3) is given by 



(2) 



»=i 



We parameterize the mixing matrix V using the Maiani representation [ ITJ | of the mixing 
matrix U in the quark sector, with CP-violating phase set to zero, i.e. 



V 



( cecp 



-s e c 7 - c 6 s 1 s 



seep \ 
cgc^ - s e s^sp s 7 C/3 



(3) 



\ SgSj — CgC^Sp — CgSj — SgC^Sfj C-yCg J 

where sg = sin 9, eg = cos 9, sp = sin/3 In Eq.(^), the angle 9 refers to mixing between 

V\ and z/ 2 , /3 to V\ and z/ 3 and 7 to z/ 2 and z/ 3 , respectively. 

In our analysis, we will discuss the limits on neutrino masses and mixing angles which 
can be inferred from purely leptonic decays of pion and leptons (/z and r). That is, the 
decay rates to be used are 
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r(vr — > eP e ), r(vr — ► ^P M ), 
— ► eP e i^), 
r(r — ► eu e u T ), r(r — ► AiP M z/ r ). 

These are the best known experimental quantities which do not involve hadrons in the 
final states, hence introducing no further unnecessary complications in the calculation of 
decay widths. We will not consider the K decays, even though some experimental determi- 
nation of its decay widths into leptons are almost as good as those of the pion. Its properties 
are quite similar to those of the pion and its data do not provide any additional (or critical) 
information. 

Here, we briefly summarize the formulas to be used. (Details are given in Appendix.) 
For the pion, the decay rate into two leptons in a general case of three massive neutrinos 
with mixing is given by 

n2f2TT2 3 3 

rfr — in) = G ; f "* m * -Kri-Y, \Vu\ 2 p- 1 , (4) 
87r i=l 

where G is the Fermi constant (see comment below), m n is the pion mass, f n is the pion decay 
constant and U U( i is the ud component of the mixing matrix in the quark sector. The matrix 
element-phase space function, P? 1 , denotes the quantity of our interest which contains part 
of the matrix element and the entire phase space, and it depends on the neutrino masses, 
rrii (i — 1, 2, 3), as well as on the pion and lepton masses. TZ n i is a factor that represents the 
radiative corrections to the process. We stress here that TZ n i depends on the pion mass as well 
as on the lepton mass. The complete expressions for P[ l and TZ wt are given in Appendix. We 
wish to emphasize here that the value of f n = (130. 7±0.1 ±0.36) MeV quoted in PDG data 
set is obtained from the decay n — > fiu^ + //iyy under the hypothesis that the neutrinos are 
massless with no mixing. For massive and mixed neutrinos, the above value represents the 
quantity [/JX)f=i ^nifPi^ / Pq^^ 2 rather than /„., where Pq^ is the matrix element-phase 
space function for massless neutrinos. For this reason, in the following, we will consider, as 
was done by PPZ, only the ratio of the two leptonic decay widths of the pion, in order to 
cancel out the dependence on the unknown quantity f n . 

The decay width for a lepton decaying into three leptons is given by 

ni' — iw) = ^| -n v ■ E (5) 

where my is the mass of the decaying particle. Again, P-'j l is the matrix element-phase space 
function which depends on the masses of all the particles involved in the decay process. The 
leading radiative corrections are denoted by IZy] they depend only on the mass of the 
decaying lepton. Also, the expressions for P\^ and 1Zi> are collected in Appendix. 

Again, it is crucial to emphasize, as noted by PPZ, that the experimental value of G M 
quoted in the PDG data set becomes the Fermi constant G only for the massless neutrinos, 
as implied by the Standard Model. Since new physics beyond the Standard Model is what we 
wish to investigate, the coupling constant G M that one measures in the muon decay should 
be interpreted as 
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where Pqq is the matrix element-phase space function for massless neutrinos. Therefore, 
in a general case of massive neutrinos with mixing, the weak coupling constant G, which 
enters in the calculation of all the weak processes like e.g. Eq(^|) and Eq.(^[), is not directly 
measured from the muon decay. The quoted number G M = (1.16639 ±0.00002) • 10~ 5 GeV~ 2 
is valid only for massless neutrinos. 

An alternative way to obtain G comes from the Standard Model of electroweak interac- 
tions, where G, the fine-structure constant a, the weak mixing angle 8\y and the W boson 
mass are related as iH|lT 



Q = ™ ft) 

V2m 2 w sm 2 e w (l - Ar)' 

The above equation relates the low-energy quantities G and a to the quantities defined 
at the electroweak scale, sin 2 9w and mw- The radiative effects are taken into account in 
Ar, which in the first order in a can be cast in the following way 

Ar = Ar Q + Ar t + Ar f + Ar 9 + Ar b (8) 

The first term in Eq. (||) describes the running of a from the zero-momentum scale to the 
electroweak scale, the second term is due to the radiative contribution of the top, Ar/ and 
Ar g are contributions coming from the leptons and quarks (except the top) and the last one 
is due to the Higgs-boson loops. The actual form of the radiative correction Ar depends on 
the renormalization scheme adopted to perform the calculations. 

Hence, in the Standard Model only three out of the four quantities G, a, sin 2 9w and 
raw are independent. Since the Fermi constant (for massless neutrinos), the fine-structure 
constant and the Weinberg angle are the best known quantities, the standard practice is to 
consider them as being independent, and then use Eq.(0) to predict the W mass. 

The aim of our analysis is to investigate the possibility that neutrinos are massive and 
mixed, therefore G, as was discussed, is different from G^ and hence should be treated 
as unknown. Namely, we reverse the standard procedure and use Eq.(^) to calculate the 
value and the allowed 1-a range for the Fermi constant, using a, sin 2 Q w and m w as input 
parameters. Obviously we expect the accuracy of the value G to be rather poor, i.e. of the 
same order of the experimental error on rriyy. It is encouraging, however, that the latest 
measurement on the W mass are at the level of 2% [O. 
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Special caution must be exercised in choosing the value of sin 2 9w to be used in Eq.(^). 
This parameter has been obtained to a very good accuracy by a number of different ex- 
periments (see, e.g. Ref. M and references quoted therein, for a summary of the different 
techniques used to extract sin 2 6^). All these methods involve the measurement of cross 
sections, and the fits of the data (from which sin 2 9w is obtained) are made using G for 
massless neutrinos (G = Gy) as an input parameter (directly in the cross section evalua- 
tions, or indirectly through the interference terms in the calculations of the asymmetries). 
This means that the best mean value of sin 2 9w quoted in the PDG data cannot be used in 
Eq.([7p, because it is not independent of G itself. Therefore, we perform our calculation of 
G in the on-shell scheme of renormalization JTT], where the following relation is defined to 
be true at all orders in the perturbation theory 

^ 2 e w = i- 7 ^-. (9) 

™>z 

In this way sin 2 9w is obtained independently of G, even though its uncertainties will 
become larger than that of the best mean value in PDG data due to the uncertainties in 
the weak boson masses. (A further comment: strictly speaking, the determination of mw 
and rriz could, in principle, be also affected by the value of G. In fact, the usual way to 
determine the weak boson masses involves a global fit of a number of physical quantities and 



G is again taken as an input parameter of the fit |L3| . Therefore, the actual value of G could 
influence very slightly the determination of mw and mz- However, the correlation between 
the value of G and the position of the poles in the W- and Z-pro duct ion cross-sections 
(basically the weak boson masses) cannot be substantial. For this reason we assume the 
quoted measured values of mw and mz as being independent of G.) 

In summary, we have chosen to determine the value of the Fermi constant and its I— a 
allowed interval by means of the Standard Model relation Eq.(|7|), in the framework of the 
on-shell renormalization scheme where sin 2 9w is defined by Eq.flj^). The actual expressions 
for the radiative corrections Ar in the on-shell scheme can be found in Ref. [|nj. The top 
contribution Ar t , not contained in the original paper by Sirlin, is || 

3a mf 1 
167r sin 2 9 W m w tan 2 9w 

Hence our input parameters are: mz, mw and a. (Obviously, the low-energy value of 
the fine structure constant, obtained from the quantum Hall effect, is independent from G). 
Some additional input parameters, which enter in the radiative corrections Ar, are the mass 
of the top quark m t and the Higgs boson mass m#. For the top mass, we use the recent 
CDF measurement m t = (176 ± 18) GeV 0, both for the 1992 and 1995 data sets. The 
Higgs mass is varied in the interval (60 GeV, 1 TeV). 

In the case of the 1992 data set, we use the 1992 PDG values for mz and mw- The 
results are (here and hereafter, the errors are propagated quadratically) 

sin 2 9 W = 0.2258 ± 0.0050 

Ar = 0.0402 ± 0.0088 (11) 
G = (1.162 ± 0.029) • 10" 5 GeV" 2 
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For the 1995 data, in order to reduce as much as possible the uncertainties in the de- 
termination of G, we use the latest data available: rriz = (91.1884 ± 0.0022) GeV [15] and 
m w = (80.410 ± 0.180) GeV 0. This gives 



sin 2 9 W = 0.2224 ± 0.0035 

Ar = 0.0396 ± 0.0088 (12) 
G = (1.174 ± 0.022) • 10~ 5 GeV" 2 

The values of G which we will use in the evaluation of the leptonic decay widths are those 
given in the two previous equations. Unfortunately, the errors on G (of the order of percent) 
are much worse than the errors on which are of the order of 10 -5 . We stress that since 
we are interested in the investigation of massive and mixed neutrinos, the method discussed 
above turns out to be the only viable way to determine the Fermi constant. Unfortunately, 
however, this procedure determines G with errors of a few percent level and thus makes the 
accuracy in the calculation of the decay widths rather poor. Nevertheless, as will be seen 
in the following, it is still possible to use them as an important constraint on the neutrino 
parameters (mainly on the neutrino masses). In fact, a deviation (if any) of G^ from the 
Standard Model value of G given in Eq.(|7]) signals the massive and mixed neutrinos. 

Measured experimentally are the branching ratios of the above decay processes. The 
branching ratios are simply related to the previously defined quantities as 

B R( 7r — lv{) = r(?r ~T Wl) = r w ■ r(7r — W t ) (13) 



B R(J' — = ~T lVlVv) = r v ■ T(l' — > Ivw) (14) 

l v 

where and Ty are the total widths of the pion and the decaying lepton, respectively, 
and r n and Ty are the corresponding lifetimes. In order to directly extract information 
about the neutrino mass and mixing, instead of the lepton decay widths themselves, we use 
the following quantities which are simply proportional to the decay widths, with common 
constants such as G and my removed, 

= a ,,BR(Z' — Ivxw) = ftp • E \Vi'i\ 2 \Vij\ 2 P!j, (15) 

where 

192tt 3 . . 

G 2 m\,Ty 

From the experimental values of the BR's and the physical quantities defined in Eqs. (|1~6|). 
we will first calculate the 1-a allowed experimental ranges for the T's. These ranges will 
then be compared with the calculated values of the T's by varying the neutrino masses and 
mixing angles. This would limit the neutrino masses and mixing angles. 

We will follow the same procedure for some ratios of the T's and the BR's. Our choice 
of the ratios is 
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In the case of lepton decays, we will perform our analysis using the two ratios, 
and The use of ratios alone is indeed simpler because uncertainties in some constant 
quantities are cancelled out, but some changes in the numerator and the denominator coming 
from phase space and mixing angles may partially be compensated. Therefore one must check 
that the calculated single partial decay widths do not lie outside the experimental allowed 
ranges. This is why we will add as an additional constraint also the three leptonic decay 
widths. Summarizing, we will use the following quantities as constraint: r , T , T , R™, 

When we evaluate the T's and the i?'s from the experimental values, we propagate the 
errors quadratically. In the calculations, we use the central values of masses of the particles 
involved. All the calculations are carried out at 1—a level. Also, it is to be noted that the 
radiative correction factors IZ's are included in the quantities T's; they are cancelled out in 
K% but not in and R™. 



III. NUMERICAL INPUTS 



We have listed, in Table 1, all the latest (1995 data) experimental inputs that will be 
used in our analysis @,|16|. In order to compare our new analysis with the previous one by 



PPZ based on the 1992 PDG data, we have also listed the 1992 PDG data § in Table 1. As 
can be seen in Table 1, the entries with asterisks signify those with noticeable changes from 
the 1992 data to the 1995 data. In particular, they include all the data on r decays and 
7i — > eU e . For the sake of comparison, in Table 1 the value of both G M and G are included. 
(We note that, strictly speaking, m T quoted in Table 1 is for massless v T . However, in this 
case, the use of m T in Table 1 does not introduce any significant modification.) 

We should notice that the errors on the T's are quite large (of the order of a few percent) 
due to the uncertainties in the Fermi constant G, as discussed in the previous Section. 
Nevertheless, the inclusion of the T constraint turns out to be effective in limiting the 
allowed intervals for the neutrino parameters, in particular the masses, as will be shown in 
the following Sections. 

Listed below are the calculated values of T's and i?'s using the latest mass values for the 
decaying particles and charged leptons under the assumption that neutrinos are massless 
with no mixing. In the case of T's, the first numbers on the right-hand side are the values 
without the radiative corrections and the second numbers represent the radiative corrections. 
All the values are in agreement with the 1995 data in Table 1 within l-a, implying that the 
1995 data are consistent with the lepton sector with massless neutrinos with no mixing. 



Calculated values for massless neutrinos with no mixing 
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(1995 data set) 




0.999813 • 0.995797 = 0.995611 
1.00-0.996 = 0.996 
0.972 • 0.996 = 0.968 



R™ = 1.233- 10~ 4 
B£ = 0.9995 
Rlf, = 1.028 



(The uncertainties in the above numbers due to the experimental errors of the quantities 
which enter in their calculations are always less than 0.01%). 

The conclusion based on the above numbers that the experimental data (1995) are con- 
sistent with the assumption of massless neutrinos with no mixing is in sharp contrast to the 
result coming from the 1992 PDG data, as obtained by PPZ. In order to further examine 
1992 data set, we have repeated the calculations using the 1992 PDG data both with and 
without the inclusion of the radiative corrections. We list, in the following, the results of 
the calculated values based on the 1992 PDG data in Table 1 with the assumption that 
neutrinos are massless with no mixing. 

Calculated values for massless neutrinos with no mixing 
(1992 PDG data set) 



= 0.999813 • 0.995797 = 0.995611 
r Te = 1.00 • 0.996 = 0.996 
= 0.973 • 0.996 = 0.969 

R™ = 1.233 • 10" 4 (R™ = 1.283 • 10~ 4 without RC) 
R!£ = 0.9995 (R!£ = 0.9998 without RC) 

K% = 1.028 (R™ = 1.028 without RC) 

(Again, the uncertainties in the above numbers are also less than 0.01%). 

Note that the calculated values of the following quantities do not lie in their corresponding 
experimental 1-a ranges: R^ e and R%%. This implies that the 1992 PDG data are indeed 
incompatible, at least within 1-a, with the assumption that neutrinos are massless with no 
mixing. If the radiative corrections are not included, Y T& is not compatible, also. 



In this Section we will first reexamine in detail the result of the PPZ analysis and then 
present the results of our new analysis. As mentioned already, the PPZ analysis was based 
on the 1992 PDG data, improved by the latest (at that time) determination of the r mass. 
In order to cancel out the dependence of the decay widths on some parameters (the pion 
mass, the quark mixing angle U u di the pion decay constant, the Fermi constant and the 
muon or tau mass), PPZ considered only the ratios of the partial decay widths; R^, Rfl 
and K!f . They did not take into account radiative corrections for these processes under the 
assumption that the radiative corrections of order of several percents are of no importance 



IV. REEXAMINATION OF PREVIOUS ANALYSIS 
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(note that the radiative corrections do not cancel each other in the ratios R* e and R^f). 
The PPZ analysis was performed in the case of one massive neutrino (z/ 3 ) and two almost 
degenerate very light (mi ~ m 2 <C m 3 ) ones. Their main result is that z/ 3 could have a mass 
in the interval 

155 MeV < m 3 < 800 MeV. (19) 

That is, all the 1992 PDG data could be fitted with m 3 in the above range. They then 
improved the upper limit in Eq.(|l9D by taking into account the constraint coming from the 
Z invisible width. The resulting allowed interval was 

155 MeV < m 3 < 225 MeV. (20) 

Given the values of m 3 inside this range, the PPZ analysis also showed that one of the 
mixing angles (namely (3) was constrained to a finite range which did not include (3 = 0. 
They obtained, for m 3 = 165 GeV, 

11.54° < (3 < 12.82°. (21) 

Although the above mass and mixing angle intervals are allowed by the ratios R, one 
must make sure that the same allowed ranges do not violate the experimental partial decay 
widths. This turns out to be the case for the mass range m 3 > 215 MeV, as will be shown 
in Section IV-A. 

In the following we will examine what would happen to the above PPZ conclusions if 
1) radiative corrections are taken into account, 2) the constraint from T's (decay widths) 
is imposed and 3) 1995 data are used together with the T constraint and the radiative 
corrections. Specifically, we will show that due to the accuracy of the present data, it 
is important to include radiative corrections and that the use of ratios i?'s alone without 
checking the partial decay widths (i.e. T's) could lead to overestimates of the allowed interval 
of the neutrino parameters. 

A. Allowed Range for Mass 

In Fig.l we have plotted the ranges for the values of m 3 which are forbidden (denoted 
by solid lines) by the ratios i?'s, by the T's and by the combination of the two. The heavy 
solid line represents the allowed region. In this plot, it is assumed as in the case of the PPZ 
that the other mass parameters are very small (mi ~ rri2 <C m 3 ) and the mixing angles 
are varied over the maximum interval (0,7r/2). The figure refers to three cases: 1992 PDG 
data, without radiative corrections (RC), denoted by 92, 1992 PDG data with RC, denoted 
by 92RC and 1995 data with RC, denoted by 95RC, respectively. For each case we present 
the three results, one with the R's, one with the T's and one with i?'s and T's combined. 
For 92, there is an allowed region for m 3 which does not include m 3 = when i?'s alone are 
used. This is the PPZ result. Most values of m 3 inside this allowed region, however, violate 
the limits on the T's, as can be seen in Fig.l. It has to be emphasized that the range for m 3 
allowed by the T's alone does not include m 3 = 0, also. 
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If we combine the two results, we find the following allowed range for m 3 

178 MeV < m 3 < 215 MeV. (22) 

The inclusion of radiative corrections changes the picture dramatically. As can be seen 
in Fig.l (92RC), the region allowed by .R's is considerably enlarged. Also the T constraint 
is modified in an important way, i.e. all the region of lower masses is now allowed, including 
m 3 = 0. When we combine the two constraints (R + T) the PDG 92 data set is consistent 
with the following finite (not including zero) mass range for m 3 

140 MeV < m 3 < 210 MeV. (23) 

signaling new physics beyond the Standard Model. It is, therefore, extremely interesting to 
repeat the analysis with the 1995 data. It is to be noted that an agreement of the above 
result with the PPZ's (Eq.(pH|)) is purely accidental. 

As can be seen in 95RC in Fig 1, i?'s alone allow two different (disconnected) regions, 
one of which includes m 3 = 0. The inclusion of r's restricts these intervals to the region of 
lower masses. The allowed regions resulting from the combined (both .R's and r's) analyses 
based on the 1995 data with radiative corrections is (for m\ ~ m 2 -C m 3 ) 

m 3 < 70 MeV and 140 MeV < m 3 < 149 MeV. (24) 

The above result is rather insensitive to the choice of m\ and rri2. For example, for 
mi < 20 KeV and m<i < 1 MeV, the above result remains unchanged. Only if rri2 is of 
the order of a few MeV all the region between m 3 = and m 3 ~ 149 MeV is allowed. 
Therefore, we conclude, based on the 1995 data, that the 1992 PDG data were internally 
inconsistent because of the poor data on r. Moreover, the accuracy of the current data or 
even the 1992 PDG data warrants the inclusion of the radiative corrections for any precision 
analysis. Furthermore, one can see that the use of the .R's alone without the constraint from 
r's can give rise to overestimated allowed regions. 

B. Allowed Range for Mixing Angles 

So far we have reexamined and discussed the PPZ using the 1992 PDG data. We have also 
carried out a similar analysis using the 1995 data with a conclusion that there is an allowed 
window for m 3 including m 3 = 0, as well as an isolated range of m 3 (for mi ~ m 2 -C m 3 ). 
In this Sub-Section we carry out a similar analysis for mixing angles. The PPZ concludes, 
based on the ratios alone from 1992 PDG data, that the angle f3 has a finite allowed range, 
11.54° < (3 < 12.82°, whereas the angle 7 is restricted to 7 < 4.05°, including zero. This 
allowed region is shown in the /3—j plane as an area filled with circles in Fig. 2a. In addition, 
the constraint imposed by the r's alone is indicated by the dotted region. (It is to be pointed 
out that in Fig.2 the mass parameters are mi ~ m 2 ~ and m 3 = 200 MeV). The region 
allowed by the two constraints combined is denoted by the dark area. 

Now, the inclusion of radiative corrections to the PPZ analysis leads to changes in the 
allowed regions of Fig. 2a. The region allowed by -R's alone is enlarged and moved towards 
the origin, but still it does not include the origin that corresponds to the case of neutrinos 
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without mixing. This is shown by circles in Fig. 2b. Also shown in Fig.2b is the allowed 
region based on T's alone, with the radiative corrections included (dots). The darker area is 
the region allowed by the i?'s and T's combined. Therefore, the common allowed region is 



7<7.2 




(25) 
(26) 



When we use the 1995 data (with radiative corrections), the situation again changes 
dramatically. Fig. 2c shows the allowed region in the j3—j plane which is obtained by using 
the ratios R alone (circles). The allowed region is shown to move farther away from the 
origin. The addition of the constraint coming from the T's completely washes out the 
region, i.e. no allowed region. This is self-evident because the mass m 3 = 200 MeV is not 
allowed by the 1995 data, as can be seen in Eq. (|24]) . Thus, even in the case of the 1995 
data, neglecting the T constraint could lead to erroneous conclusions. 



In this Section we present the results of a more detailed analysis of the limits that can be 
set on the neutrino masses and mixing angles by using the 1995 data and with the inclusion 
of radiative corrections. It is to be pointed out here that a complete, combined analysis of 
the masses (mi, m 2 and 777,3) and mixing angles (8, (3 and 7) is very much involved and is 
beyond the scope of this paper. Even the presentation of results of such an analysis would 
be problematic. Therefore, we have simplified the analysis by fixing some parameters and 
varying others. In order to see various correlations among the masses and the mixing angles, 
we present several allowed regions in the two dimensional plots for several combinations of 
mass and mixing angle. 

First, we present the absolute upper limits on the three neutrino masses, mi, m 2 , and 
m 3 , independent of the values of mixing angles, 9, (3 and 7. We have obtained these limits 
by varying the mixing angles over the entire interval between and 7r/2 and by taking into 
account all the constraint which we have discussed in the previous Section (i?'s and T's) 
and the radiative corrections. (It should be stated that we are not carrying out a statistical 
analysis of all the relevant data.) Instead of three dimensional plots, we present, in Fig.3, 
the allowed region (dotted area) in the mi-m 2 plane and the allowed region in m 2 -m 3 plot 
in Fig. 4. From these plots, we can set the following absolute upper limits on the neutrino 
masses (based on Y and R constraint, at 1-er level) 



where the limits on mi and m 2 are mainly due to the i?'s, whereas the limit on m 3 comes 
from the V constraint. The limit on m 3 has already been mentioned in the previous Section. 
Although the limit on m x is rather poor, the limit on m 2 is larger by a factor of less than 
thirty than the latest limit from the kinematical analysis of the ir — > fi + decay. Similarly, 
the limit on m 3 is larger only by a factor of six. It is quite interesting that the accuracy 



V. ANALYSIS AND RESULTS 



mi < 100 KeV, 
m 2 < 7.5 MeV, 
m 3 < 149 MeV. 



(27) 
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of the present data on the decay rates and branching ratios is already sufficiently good 
enough to set limits on m 2 and m 3 which agree, in one order of magnitude, with the results 
from more involved kinematical determination. The important difference between the upper 
limits given in Eq . (|27|) and those in Eq.(|I]) is that the former is valid independently of 
mixing angles whereas the latter is valid only for the case of no mixing. Furthermore, the 
improvement in the data from 1992 is obvious from the conclusion that the upper limit on 
m 3 is set to 149 MeV and the limits are consistent with massless neutrinos with no mixing, 
implying the internal consistency of the data. 

Next, we discuss some correlations among the masses and the mixing angles. The first 
example to be presented is the allowed region in the m 2 -sin 2 6 plot shown in Fig.5. In this 
plot, we have set m 3 to be 24 MeV and u 3 is assumed to be decoupled, due to its heavy 
mass, from v\ and z/ 2 so that (3 = 7 = 0. Also, for definiteness, we have taken mi = 5 
eV, but the conclusion remains unchanged as long as mi is less than ~ 20KeV. The solid 
and dashed lines delimit allowed regions based on the use of T's and the constraint from 
i?'s, respectively. The allowed area is denoted by dots. Fig.5 shows that the low-angle 
regime is constrained mainly by whereas R™ is more effective in limiting the large- 
angle area (say sin 2 6 > 10~ 4 ). The allowed region shown in Fig.5 is insensitive to values 
of mi and m 3 , as long as they are m x < 20 KeV and m 3 < 50 MeV. The so-called small- 
and large- angle solutions of the solar neutrino deficit based on the MSW effect, which 
requires m 2 — m\ ~ 6 x 10~ 6 eV 2 and sin 2 (26) ~ 7 x 10~ 3 , and m 2 , — m\ ~ 8 x 10~ 6 eV 2 and 
sin 2 (2#) ~ 0.6, respectively, are well within the allowed region in this plot. 

In the next example, we assume that v\ is too light to couple with z/ 2 and z/ 3 . That is, 
only z/ 2 and z/ 3 are mixed with angle 7. We have also set m\ <C m 2 = 270 KeV. Fig. 6 shows 
the allowed region in the m 3 -sin 2 7 plane. Here, the small- and large- angle regions are 
constrained by R!j%, whereas the intermediate region (sin 2 7 ~ 10 _4 -10 -2 ) is constrained by 
R%p. The neutrino oscillation solution of the atmospheric neutrino anomaly as observed by 
Kamiokande and others, which favors m\ — m\ ~ 10" 2 eV 2 and sin (26) ~ 1, is well within 
the allowed region in Fig. 6. The size and the shape of the allowed region are insensitive to 
the assumed values of mi and m 2 as long as mi < 20 KeV and m 2 < 1 MeV. 

Although the case of a V\- u 3 mixing is unnatural in the framework of natural hierarchy 
of neutrino masses, we present it in Fig. 7. In this Figure the allowed region in the m 3 - 
sin 2 /? plane is shown. Here, as in the previous case we have set mi <C m 2 = 270 KeV. In 
this case, the most restricted limit is always imposed by R^. Again, the allowed region 
does not change significantly as long as m 1 < 20 KeV and m 2 < 1 MeV. 

In the previous figures Figs. 5-7 the limits were obtained in the special cases in which 
two of the three mixing angles were kept fixed at zero, i.e. only one pair of neutrinos is 
mixed. To show how sensitive these limits are to the fixed angles and correlations among 
the limits, we present in Figs. 8-10 the cross sections (with one angle fixed) of the three-angle 
parameter space for fixed values of the masses, i.e. sin 2 7- sin 2 /? in Fig.8, sin 2 7-sin 2 # in 
Fig. 9 and sin 2 /?- sin 2 6 in Fig. 10, respectively. In each figure, the values of the masses are 
fixed as mi -C m 2 = 270 KeV and m 3 = 24 MeV. It is interesting to note that for the above 
set of the masses, the most severely constrained angle is /?, whereas the least constrained is 
7. The limits are: 
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sin 2 9 < 3.6 • 1(T 3 

sin 2 /3<4.6-10- 7 (28) 
sin 2 7 < 7.0- 10- 3 . 

As we decrease the values of the masses used, the allowed region in each figure in- 
creases, eventually covering the entire space. Hence, no meaningful limit can be obtained, 
as expected. In order to demonstrate this sensitivity, we have shown in Figs. 8-10 the ex- 
tended allowed regions (bounded by the dashed lines) for the following values of the masses: 
mi <C m 2 = 10 KeV and m 3 = 1 MeV. In this case, the limits on the mixing angles 
substantially increase to 

sin 2 9 < 1 

sin 2 (3 < 2.6 • 10~ 4 (29) 
sin 2 7 < 1. 



VI. SUMMARY AND CONCLUSIONS 

We have analyzed and compared the 1992 and 1995 data on the 7r, fi and r decays in the 
framework of three generations of massive neutrinos with associated mixing. First we have 
confirmed the surprising result of PPZ based on the 1992 PDG data that when only the 
ratios, i?'s, are used without radiative corrections, the 1992 data are inconsistent with the 
picture of massless neutrinos with no mixing, signaling new physics beyond the Standard 
Model! More specifically, PPZ has shown that when m x and m 2 are assumed to be much 
less than m 3 , m 3 is found to be within the interval 155 MeV < m 3 < 225 MeV and the 
v\ - v 3 mixing angle, (3, in the interval, 11.54° < (3 < 12.84°. This isolated allowed region 
survives even if we introduce the T constraint. 

Again, the 1992 PDG data set with radiative corrections reproduces the results that agree 
qualitatively with those of PPZ. The allowed range of m 3 is 140 MeV < m 3 < 210 MeV and 
the allowed mixing angles are 4° < (3 < 10° and 7 < 7.2°. This clearly shows that the 1992 
PDG data set suggests massive and mixed neutrino. 

In order to see if this rather surprising result still remains valid or not with the improved 
data of 1995, we have carried out a comprehensive analysis of the 1995 data by using both 
the ratios, i?'s, and decay widths, T's, and by including the radiative corrections. The 1995 
data are shown to be internally consistent, and furthermore consistent with the picture of 
massless neutrinos with no mixing. 

Limits on the masses derived from the analysis are as follows: 

m 1 < 100 KeV, 

m 2 < 7.5 MeV, (30) 
m 3 < 149 MeV. 

These bounds on the masses are such that the imposed constraint (-R's and T's) are 
fulfilled in their 1-a intervals. 
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Although the above limits are less stringent than those in Eq.(|l]) from kinematical de- 
terminations, it is important to note that the limits given in Eq.(30) are completely in- 
dependent of mixing angles. Therefore, if neutrinos are massive and mixed, their masses 
can be heavier than the limits given in Eq.(|I|). Of course, such heavy neutrinos, if sta- 
ble, are not allowed by the well-known cosmological limits, Z)i^( z/ i) < 20 ~ 30 eV (with 
h = (Hq s Mpc)/100 Km ~ 0.5, where Ho is the Hubble constant). If unstable, the decay 
of neutrinos must be such that it should not disturb the standard nucleosynthesis scenario 
and should not violate the observed limits on cosmic electromagnetic wave backgrounds. 

It is not possible to obtain absolute limits on the mixing angles, because they strongly 
depend on the input values of the masses. Examples of limits on the angles for definite 
values of the masses have been derived and reported in the previous Section. 

To conclude, it is gratifying that the accuracy of the current data is already good enough 
to set limits on m 2 and m 3 which agree, in the order of magnitude, with the results from more 
involved kinematical determination, although the limits on the mixing angles are still rather 
poor. Further improvements of the data on r, my/ and mz in the future may significantly 
improve the limits. 
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APPENDIX A: FORMULAS FOR PION AND LEPTON DECAYS 

In this Appendix we present, for completeness, the expressions for the pion and lepton 
decays for the general case of three massive neutrinos with mixing. Radiative corrections to 
these processes are included with a brief comment. 

1. Pion decay 

The partial decay rate of the pion into two leptons is given by 

Q2 f2jj2 3 3 

r(vr — > w t ) = u J ^ m - • n wl ■ £ \v H \ 2 P[ l , (Ai) 

where the phase space-matrix element factor, Pf l , is given by 

Pf = 6(m n - m, - raj ■ K + C - K - <&) 2 ] A 1 ^!, <£, €)• (A2) 



15 



In the above, G is the Fermi constant, is the pion decay constant, U is the mixing matrix 
in the quark sector, 



Sin 



mi 



1,2,3) 



(A3) 
(A4) 

(A5) 

The quantity TZ^i in Eq. ( |A1|) describes the leading radiative corrections to the pion decay 
process (|,[17| given by 



and A is the standard kinematical function 

.2,2, 2 



X(x, y, z) — x + y + z — 2{xy + yz + xz). 



2a 



m z 



1 + — In — 



7T 



X 1 



a 

TV 



m r 



In 



m. 



1 + -F(5 l7T ) 

7T 



+ Ci + cAln 

nip 



m r 



mr 



3 — 2 



+ 



(A6) 



where 



F(x) = 3 In x 
-2( 



13 - 19x 2 8 - 5x 2 



8(1 -x 2 
1 + x 2 



2(l-x 



2\2 



x lux 



x. 



\nx + 1 ln(l -x 2 ) + 2 



1 + x 2 
1-x 2 



L{l-x 2 ) 



(A7) 



Here, m p = 796 MeV is the p meson mass, mz the Z boson mass and a is the fine structure 
constant. Also, in the above, L(z) is defined by 



Liz) 



ln(l - t) 



dt. 



(A8) 



The first square bracket in Eq. (|A6|) represents the electromagnetic short-distance cor- 
rection. Its value is slightly modified when higher-order effects and QCD corrections are 
taken into account, i.e. |T?| 



1 H In 

7T 



m z 



1.0232. 



(A9) 



The second and the third brackets denote the QED corrections to the decay of a point- 
like pion |T7|JT8|| . Following a general practice, we neglect the terms with the CVs, whose 
numerical values have large uncertainties [17]. Therefore, we use the following simplified 
expression 

3a 



K*i = 1.0232 



a 



1 + -F(6 lv ) 

7T 



1 



In'™' 



2tt \m 

Note that TZ n i depends both on the pion and the lepton mass. 



(A10) 
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2. Lepton decay 



The partial decay rate of a lepton into three leptons is given by 



G 2 mf, ' 

1927T 3 



T(l' — W M ,) = ^ -Kv-Y, \Vn?\V^I%\ (All) 



where 



and 



= ( mj< _ mi _ m . _ m .) . 2 / (a; 2 - 4<Jg / ) 1 / 2 - Mdx (A12) 



= ^ 1/2 [x(l + 4 - x)A + (2-x)(x- 261)8], (A13) 
1 



A 



B 



(1 + 81 - xf 

x [(1 + 8 2 , - xf - 2(1 + 81, - x)(4 + 8%,) + (81, - 5%,)\ (A14) 
1 

(1 + 52, - x) 2 

x [(1 + 81 - xf + (1 + 81 - x){5l, + 5%) - 2(4 - 4,) 2 ]- (A15) 

In the above 

8 IV = ^ (A16) 
m v 

Tfl ' 

8 iV = — (i = l,2,3) (A17) 

X min = 28 lV (A18) 

x max = 1 + 8 2 w - (8 iV + 5 jV f. (A19) 

The quantity IZy describes the leading radiative corrections to the lepton decay process 
[ 19| which are given by 



a /25 

H 71' 

2ir V 4 



i + iH- p (A20) 



where m\y is the boson mass. Note that Hp depends only on the mass of the decaying 
lepton, but not on the mass of decay products. 
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FIGURES 



FIG.l Allowed (thick solid lines) and forbidden (solid lines) intervals of m 3 by -R's, r's 
and .R's and r's combined, respectively, for 1992 and 1995 data sets (mi ~ m 2 <C m 3 is 
assumed). Mixing angles are varied over the interval (0, 7r/2). Label 92 refers to the result 
of 1992 PDG data, without radiative corrections in the calculation; 92RC and 95RC denote 
the results of 1992 PDG and 1995 data sets with radiative corrections, respectively. The 
mass scale is in arbitrary units. 

FIG. 2 Allowed (dotted) regions in the f3-^ plane for 1992 and 1995 data sets with 
m\ ~ m 2 — and m 3 = 200 MeV: (a) allowed region by .R's alone (circles), r's alone (dots) 
and R + T combined (dark area) of 1992 PDG data, (b) allowed region by -R's alone (circles) 
and by r's alone (dots) of 1992 PDG data with radiative corrections (the darker area is the 
region allowed both by the r's and the -R's), (c) allowed region by .R's alone (circles) of 1995 
data with radiative corrections. 

FIG. 3 Allowed (dotted) region in the m r m 2 plane for 1995 data with radiative correc- 
tions. Both .R's and T's constraints are imposed. Mixing angles are varied over the interval 

(0.7T/2). 

FIG. 4 Allowed (dotted) region in the m 2 -m 3 plane for 1995 data with radiative correc- 
tions. Both i?'s and T's constraints are imposed. Mixing angles are varied over the interval 

(0.7T/2). 

FIG. 5 Allowed (dotted) region in the m 2 -sin 2 9 plane for 1995 data with radiative 
corrections and mi = 5 eV, m 3 = 24 MeV and f3 — 7 = 0. The regions below the solid and 
the dashed lines are allowed by T's alone and .R's alone, respectively. 

FIG. 6 Allowed (dotted) region in the m 3 -sin 2 7 plane for 1995 data with radiative 
corrections and m x = 5 eV, m 2 = 270 KeV and 6 — (3 — 0. The regions below the solid and 
the dashed lines are allowed by T's alone and .R's alone, respectively. 

FIG. 7 Allowed (dotted) region in the m 3 -sin 2 /3 plane for 1995 data with radiative 
corrections, and mi = 5 eV, m 2 = 270 KeV and 6 = 7 = 0. The regions below the solid and 
the dashed lines are allowed by T's alone and -R's alone, respectively. 

FIG. 8 Allowed regions in the sin 2 7-sin 2 f3 plane for 1995 data with radiative corrections. 
The dotted area inside the solid line is for mi <C m 2 = 270 KeV and m 3 = 24 MeV. The 
area on the left of the dashed line is for m x <C m 2 = 10 KeV and m 3 = 1 MeV. Note the 
enlargement of the allowed region as masses decrease. 

FIG. 9 Allowed regions in the sin 2 7-sin 2 6 plane for 1995 data with radiative corrections. 
Notation and interpretation are the same as in Fig. 8. 

FIG. 10 Allowed regions in the sin 2 /3-sin 2 6 plane for 1995 data with radiative corrections. 
Notation and interpretation are the same as in Fig. 8. 
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TABLES 

TABLE I. List of 1992 and 1995 data sets. The entries with asterisks denote data which have 
been improved by a significant amount. 







1992 data 


1995 data 






(0.51099906 ± 0.00000015)MeV 


same as 1992 




"V 


(105.658389 ± 0.000034)MeV 


same as 1992 




m T 


(1784.1 ± 3.6)MeV 


(1776.96 ± 0.31)JVleV 




777^ 


(139.5679 ± 0.0007)MeV 


(139.56995 ± 0.00035)MeV 




T M 


(2.19703 ± 0.00004) • 10 b s 


same as 1992 






/nor - i r*\ ~\ r\ — 1 5 

(305 ± 6) ■ 10 s 


(291.6 ± 1.5) -10 s 




777^ 


(91.173 ± 0.020) • 10 d MeV 


(91.1884 ± 0.0022) • 10 d MeV 




771 W 


(80.22 ± 0.26) • 10 MeV 


/on din i r\ i on\ i n3 a it \t 

(80.410 ± 0.180) • 11) MeV 






(1.16639 ± 0.00002) • 10 ii MeV 


same as 1992 




G 


(1.162 ± 0.029) • 10~ n McV~ 2 


(1.174 ± 0.022) • 10~ n McV~ 2 


St 


iiK(7r — ► ez/ e J 


(1.218 ± 0.014) • 10 


(1.230 ± 0.004) • 10 4 




BR(7T ► /T.7^) 


/~\ riAAiTTrin i i"i r\r\r\r\r\-\ A 

0.9998782 ± 0.0000014 


0.9998770 ± 0.0000004 




BR(/x — ► ez/ e i/ M ) 


1 


same as 1992 


* 


1)1) /' , _ — , , \ 

dK(t — ► ev e v T ) 


0.1/93 ± U.U026 


0.1 r /9 ± U.U0U9 




t>t~) /" — ^ 
BR(r — ► hvuVt) 


0.1758 ± 0.0027 


0.1733 ± 0.0009 


* 




1 003 ± 050 


983 ± 037 


* 


r e 


0.944 ± 0.053 


0.979 ± 0.037 


* 




0.925 ± 0.052 


0.954 ± 0.036 


* 




(1.218 ±0.013) • 10~ 4 


(1.230 ±0.004) • 10~ 4 


* 




1.063 ± 0.028 


1.0038 ± 0.0073 


* 


pre 


1.020 ± 0.022 


1.0265 ± 0.0074 
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